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Charge dynamics in a dimer Mott insulating system, where a non-polar dimer-Mott (DM) phase and 
a polar charge-ordered (CO) phase compete with each other, are studied. In particular, collective charge 
excitations are analyzed in the three different models where the internal-degree of freedom in a dimer is 
taken into account. Collective charge excitation exists both in the non-polar DM phase and the polar CO 
phase, and softens in the phase boundary. This mode is observable by the optical conductivity spectra 
where the light polarization is parallel to the electric polarization in the polar CO phase. Connections 
between the present theory and the recent experimental results in /t-(BEDT-TTF)2Cu2(CN)3 are discussed. 
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Exotic ferroelectricity beyond the conventional displacive- 
type and order-disorder-type classifications has been recog- 
nized as one of the central issues in recent condensed mat- 
ter physics. Multiferroics, i.e. coexistence of ferroelectric- 
ity, ferroelasticity and ferromagnetism, are one of the ex- 
amples. Non-collinear spin orders in Mott insulators break 
the space-inversion symmetry through the symmetric- and 
antisymmetric-exchange interactions. In this sense, the se- 
ries of materials is classified as a spin-driven ferroelectricity. 
Another example is seen in the so-called quarter-filled sys- 
tem. At low temperature, an electronic charge order occurs 
so as to break the space-inversion symmetry. This is termed 
charge-driven ferroelectricity or electronic ferroelectricity. ,2) 
Several transition-metal oxides and organic compounds have 
been studied intensively and extensively as candidates of the 
electronic ferroelectricity. Characteristics in electronic ferro- 
electricity are dielectric fluctuation. Since the electron mass 
is much smaller than the ion mass and the proton mass, large 
dielectric fluctuation due to electronic charge dynamics is ex- 
pected to play a crucial role on the ferroelectricity. 

A layered organic compound k-(BEDT-TTF)2Cu2(CN)3 
and related materials are recognized as candidates of elec- 
tronic ferroelectric materials. A crystal lattice consists of 
cation and anion layers stacked along the z axis. In a cation 
layer, two BEDT-TTF (ET) molecules form a dimer where 
bonding and antibonding orbitals are constructed from the ET 
molecular orbitals (MO). Since one hole exists per dimer in 
the chemical formula, the antibonding MO is occupied by 
one hole. Thus, this system is identified as a Mott insula- 
tor, termed a dimer-Mott (DM) insulator, in the strong dimer- 
izaion limit. 3 ' 4) No long-range spin order down to few mK 
suggested experimentally has attracted a number of researchs 
in the viewpoint of realization of spin-liquid state. 5 ~ 7) Re- 
cently, an anomalous enhancement in the dielectric constant is 
discovered around 30K. 8) This observation triggers reconsid- 
eration of the DM insulator picture and suggests a possibility 
of the polar charge distribution inside dimers and electronic 
ferroelectricity. 9 ~ 12) 

In this Letter, charge dynamics in a DM insulating system 
are examined. In particular, the collective excitations in the 
non-polar DM phase and the polar charge ordered (CO) phase 
are focused on. Collective excitation is the central concept 



not only in ferroelectricity but also in other symmetry broken 
state. 13 ~ 15) The soft-phonon mode in the displacive-type fer- 
roelectricity, the soft-proton mode in the hydrogen-bond type 
ferroelecticity, and the magnon induced by electric field in 
multiferroics are typical examples. We propose that the col- 
lective charge excitation also exists in the DM insulating sys- 
tem, where the non-polar DM phase and the polar CO phase 
compete with each other. This mode represents the polariza- 
tion flip in the polar CO phase and the intra-dimer bonding- 
antibonding excitation in the DM phase, and softens in the 
phase boundary. This soft collective mode is observable by 
the optical conductivity spectra where the light polarization 
is parallel to the electric polarization in the polar CO phase. 
Connection between the present theory and the recent experi- 
mental results in k-(ET)2Cu2(CN)3 16 ' 17) are discussed. 

To examine the charge dynamics without ambiguity, we an- 
alyze the three different models, the extended-Hubbard model 
(EHM), the spin-less Vt model (VtM) and the pseudo-spin 
model (PSM), where two molecules denoted by a and b are 
introduced in each dimer located in a crystal lattice. The EHM 
is defined as 
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where c, MS is an annihilation operator for a hole at the i- 
th dimer with spin s(=^,|) and molecule /x(= a, b), and 



s ) is a number operator. We 



introduce the intra-molecule Coulomb interaction ([/), the 
intra-dimer transfer (£a), the intra-dimer Coulomb interac- 
tion (Va), the inter-dimer transfer (t^j 1 ), and the inter-dimer 

Coulomb interactions (V^ ). When we focus on the charge 
dynamics where the excitation energy is less than the effec- 
tive intra-dimer Coulomb interaction U e s = [U + Va+ 4i a — 
y/(U -V A ) 2 + 16t A }/2, the spin-less VtM is valuable as an 
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effective model. This is defined by 
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where /i M is an annihilation operator for a spin-less fermion 
at the i-th dimer with molecule /i, and n' ifl (= /jL/i/j) is a 
number operator. The transfers and the Coulomb interactions 
for fermion in the VtM are defined in the same way with those 
in the EHM. 

When the fermion number in each dimer is fixed to be one, 
the charge state in a dimer is described by the pseudo-spin 
(PS) operator with an amplitude of 1/2. This is defined as 

Qi = (V 2 ) £ MM ' fi^n^'fi^ 1 with the Pauli matrices a and 
h = £* U^fa where U = (l/v^X^ + a x ). The eigen 
states for Qf with the eigen values of 1/2 and —1/2 are the 
charge polarized states \a) = //jO) and \b) = f} b \0), respec- 
tively, and those for Q\ with 1/2 and —1/2 are the bond- 
ing state = (| a) + |6))/v2 and the antibonding state 
| a) = (| a) — \b))/y/2, respectively. When we assume in the 
VtM that the inter-dimer transfers and Coulomb interactions 
are smaller than the intra-dimer Coulomb interactions, we ob- 
tain the low-energy effective Hamiltonian up to the orders of 
OCH t 2 )andOCH{,)as 
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where the first and second terms originate from the first and 
last terms in Eq. (2), respectively, and Wij (= V^ a + V-j — 
— Vfj a ) is the effective inter-dimer Coulomb interaction. 
The third term represents the exchange interactions, which are 
given by 
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where we introduce Qf = ± iQf and mf = ^±Q*. The 
exchange constants are given in Ref. 18) We note that a related 
model in a one-dimensional system is analyzed in Ref. 19) 

The EHM and the VtM are analyzed by the combined 
method of the exact-diagonalization based on the Lanc- 
zos algorithm, and the mean-field approximation. We in- 
troduce the mean-fields, which act on the edge-sites of 
the cluster, by decoupling the inter-dimer Coulomb inter- 
actions, such as riifiTijft' -> (nif^rijfti + n ia (ni 
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The mean fields are determined to be consis- 
tent with the electronic state inside the cluster. The charge 
dynamics are examined by calculating the optical absorp- 
tion spectra defined by a^(u) = — (e 2 /N)lm(0\j^(u — 
H + Eq + irfl^j^O), and the PS dynamical correlation 

functions defined by iV c (±) (q,w) = -lm(0|i^ (±) (a; - 

H + Eq + ir?) _1 i_q |0), where £ takes a two-dimensional 
Cartesian coordinate, |0) and Eq are the ground-state wave- 
function and energy, respectively, r\ is an infinite decimal 
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Fig. 1. (a) A schematic lattice structure for k-(ET)2Cu2(CN)3. Ellipses 
represent the ET molecules, and two molecules in a dimer are indicated 
by a and b. Black and gray lines represent the transfers and the Coulomb 
interactions introduced in the EHM and the VtM, and symbols A, B, p and 
q correspond to the subscripts in t and V in the Hamiltonians (see text). 
Two inequivalent dimers in a unit-cell are denoted by Dl and D2. We 
take lattice constants a x = 1 and a y = \/3, a distance between the two 
ET molecules in a dimer d = 0.5, and the angle of the molecular plane 
<j> = 7r/4, Shaded circles and a thick arrow represent schematic hole dis- 
tributions and a direction of the electric polarization in the polar CO phase, 
respectively, (b) The DM order parameter — (Q5_q) (soild lines) and the 
CO parameter (Q^_ ) (dash lines) for the 8-dimer cluster of the EHM, 
and those for the 12-dimer cluster in the VtM. Parameter values are chosen 
to be U = 8, V A = 2 for the EHM and V A = 2 for the VtM. 



constant, and N is the number of dimers in a cluster. We 
define the current operators j = iX}(ijW' 



,) for the EHM and j 



) (/iUv - 4' f*J for the VtM ' 
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and 4 (±) = (1/2V) (E ie Di±E. ie D 2 ) QlcM-iq ■ Ri) 
where X^6_di(D2) re P resen ts a summation for the site i which 
belongs to the sublattice D1(D2) and R l = (R la + R ib )/2. 
A schematic arrangement of molecules and interactions are 
shown in Fig. 1(a). By taking into account the relations of 
tA > *b — t p > t q in k-(ET) 2 Cu 2 (CN) 3 , we set the 
intra-dimer transfer tA as a unit of energy, and t q = and 
tp = ts(= t), for simplicity. We consider the 1/r-type de- 
pendence for the Coulomb interactions, and set V q = 0.7V P 
and V B = 0.7V P for the EHM, and V q = Q.W P and V B = 
0.8V P for the VtM as a unit of V p {= V). 

First we show the ground state properties. In Fig. 1(b), we 
show the DM order parameter — (Qg =0 ) and the polar CO 
parameter (Q^_ Q ) in the EHM. With increasing V, there is 
a critical V(= V c ) where — (Q^_ Q ) decreases and (Q^_ ) 
becomes finite from zero. The phase transition from the DM 
phase to the polar CO phase occurs at V = V c . The charge 
configuration in the polar CO phase is illustrated in Fig 1 (a) 
where the electric polarization appears along the y direction. 
Numerical results in the VtM [see Fig. 1(b)] are qualitatively 
the same with those in the EHM. 

Now we present the charge dynamics in the EHM. We be- 
gin with a simple two-dimer cluster. The optical absorption 
spectra for various V are shown in Fig. 2(a). In the DM phase, 
mainly four peaks, termed A, B, D and E, appear for small 
V. The lowest-energy peak A, shows a softening around the 
phase boundary. In the polar CO phase, the peak B splits, and 
intensities of the peak D and E are reduced. The peaks B, D, 
and E are identified as the inter-dimer charge excitations, de- 
scribed as 'Dj'Dj — > T>?T>j, where 2?" represents a state at the 
i-th dimer with n holes. On the other hand, the peak A is the 
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Fig. 2. (a) The optical absorption spectra for the EHM of a 2-dimer cluster. 
Parameter values of V are changed from 0.4 to 4 from the bottom to the 
top, and the DM and polar CO phases are realized in V < Vc = 2.39 and 
V > Vc , respectively. Other parameters are chosen to be U = 6, Va = 4, 
t = 0.5 and r\ = 0.01. (b) Schematic intra-dimer and inter-dimer charge 
excitations. Symbols a and ft represent antibonding and bonding orbitals, 
respectively, and A, B, C, D and E correspond to the peaks in (a), (c) and 
(d) (see text), (c) The optical absorption spectra a x (u) (the upper panel) 
and oiy (lu) (the lower panel) for the EHM of an 8-dimer cluster. Parameter 
values of V are changed from 1 to 2 from the bottom to the top, and other 
parameters are chosen to be U = 8, Va = 2,t = 0.6 and r] = 0.01. (d) 
The low energy parts of ctx,y (<*>) in (c). 

intra-dimer charge excitation, T>}T>j — > T>}*Dj, where 2?™* 
represents the excited state of X>". We note that the peak A is 
connected to the collective excitation of the present interest. 

The T> 2 states in a single dimer are classified as the triplet 
states, \T) = {\a t /3 t ), \a&), (|a t ^) + l«^»/v^>, 
the singlet state, \S) = (|a T /?4.) - \a^))/V2, and the 
states where two holes occupy the same MO as, \D±) = 
C±\a t ai) ± C T \PfP±). Eigen en ergies are E T = Va , E s = 
U, and E D± = [U + V A ± y/(U - V A ) 2 + W%\/2, re- 
spectively, and coefficients are given by CL/C+ = (U — 
V A )/[2E D+ - Ua -(U + Va)}. The peaks B, D, and E origi- 
nate from the inter-dimer excitations where the final T> 2 states 
are \DJ), \S) and \D+), respectively. Since the two holes oc- 
cupy the same MO and the different two MO in \ D±) and \ S), 
respectively, the peaks B and E are termed the Hubbard exci- 
tations, and the peak D is termed the inter-MO excitation [see 
Fig. 2(b)]. 

The results in a 8-dimer cluster for the EHM are pre- 
sented in Figs. 2(c) and (d). We identify the peaks A-E in 
the same way with the results in the two-dimer cluster, al- 
though multiple-peak structures are induced by the lattice ef- 
fect. Two-excitation modes, termed A + and A_, for the intra- 
dimer excitation appear. The lower- (higher-) energy mode 
is seen for the y(x) polarization, and the lower-energy mode 
shows softening around the phase boundary. 

Detailed low-energy excitations are examined in larger 
clusters of the VtM. The optical absorption spectra in a 12- 
dimer cluster are shown by red lines in Fig. 3(a). Spectra 
around 1.1 < co < 1.4 in a x (co) and 0.4 < co < 1.2 in 
Oi y (co) are attributed to the intra-dimer charge excitations, cor- 
responding to the peak A_ and A+, respectively, in Fig. 2(c) 
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Fig. 3. (Color online) (a) The optical absorption spectra a x (uj) (the upper 
panel) and a v (w) (the lower panel) for the VtM of a 1 2-dimer cluster (red 
lines). Parameter values of V are changed from 0.2 to 2 from the bottom 
to the top, and the DM and CO phases are realized in V < Vq = 1.13 
and V > Vc, respectively. Other parameters are chosen to be Va = 2, 
t = 0.8 and r\ = 0.01. The dynamical PS correlation function N x + ^ (q = 
0, lj) are also plotted by blue lines for comparison, (b) The dynamical PS 
correlation functions N x (q, ui) in the energy-momentum plane near the 
phase boundary (V = 1). Other parameter values are the same with those 
in (a). 

and (d). Remarkable polarization dependence and softening 
around the phase boundary are observed as seen in the re- 
sults for the EHM. The dynamical PS correlation functions, 
(q = 0, co) are compared with a XtV (co) in Fig. 3(a). The 
peak positions mNi~\o,co) and N x + \o, co) almost coincide 
with those in the low-energy modes of a x (co) and a y (co), re- 
spectively. Since the off-diagonal PS operators describe the 
electronic-state change inside a dimer, this coincidence im- 
plies that the lower-energy modes in a^(co) are attributed to 
the intra-dimer charge excitations. In the CO phase, the inten- 
sity of Nx (q = 0,co) decreases with increasing V, because 
the direction of PS is changed from Q z to Q x . The peaks at 
co = in N x + \q = 0,co) are the superlattice Bragg peaks 
due to the CO. The momentum dependence of N x ± \q,co) 
near the phase boundary is presented in Fig. 3(b). A clear dis- 
persion seen in the low-energy mode suggests that this mode 
is a collective charge excitation. 

To examine the low-energy excitations without the finite- 
size effect, the PSM is analyzed by the spin-wave approx- 
imation. The PS order parameters are adopted as (Qq) = 
-\5 q =o and (Q*) = for the DM phase, and (Q z q ) = 
~^6 q= ocos9 and (Q q ) = ^q^osin^ for the polar CO 
phase, where the angle 8 is determined to minimize the 
ground state energy. By applying the Holstein-Primakoff 
transformation and the Bogoliubov transformation to the PS 
operators, two kinds of bosons originating from the two- 
inequivalent dimers in a unit cell are introduced. The dis- 
persion relations for the two modes are obtained as e q ±%> = 
s/UIq ± JBq) 2 - {Jc q ± Jdq) 2 , where J Aq , Jb q , Jc q 
and Jo q are given in Ref. 20) In this scheme, the dynam- 
ical PS correlation functions are given by N^lq, co) = 
[coshr?^' — sinh0g ]S(co — e q ), where cosht?q ± ' ) = 

(J Aq ± Jsq)/4 ±) and sinhOq^ = (J Cq ± JDq)/tq • K is 
shown that N x ^ (q, co) takes peaks at e q . The e q + ^ and e q ^ 
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Fig. 4. (Color online) (a) The energy dispersions of the PS waves. Bold and 
broken lines are for e q + ^ and e g respectively. Red, green and blue lines 
show the results in the DM phase (V = 1), at the phase boundary (V = 
1.2 = V c ), and in the CO phase (V = 2), respectively. The inset shows 
the first Brillouin zone for the rectangle unit cell shown in Fig. 1(a), and 
arrows indicate the wave vector trajectories. Parameter values are chosen 
to be Va = 3 and t = 0.5. (b) The energy of the acoustic PS wave, e^t-Q, 
in the DM and CO phases. Schematic PS motions for the collective mode 
are shown in the inset. 

modes correspond to the acoustic and optical PS waves, re- 
spectively. The energy dispersions for the PS waves are shown 
in Fig. 4 (a). The acoustic mode softens around the phase 
boundary and e^o becomes zero at the boundary. 

Through the systematic analyses of the EHM, VtM and 
PSM, we obtain a whole picture for the collective charge dy- 
namics in the DM insulating system. The low-energy excita- 
tions in a x (u) and a y (ui) are identified as the acoustic and 
optical PS waves, respectively. The finite-excitation energy 
for the peak A + observed in a x (ui) at the phase boundary 
(see Fig. 2 and Fig. 3) is due to the finite size effect, and the 
energy of the acoustic mode becomes zero at the phase bound- 
ary. This soft collective mode is connected to the electric- 
polarization flip in the CO phase and the intra-dimer bonding- 
antibonding excitation, so-called the dimer-excitation, in the 
DM phase. The results are summarized in Fig. 4(b). 

Finally, the recent experimental results are interpreted from 
the view point of the present collective charge excitation. 
The optical conductivity spectra have been measured in k- 
(ET) 2 Cu 2 (CN) 3 , 17 - 21) and other K-type ET compounds, 22 24 ) 
and studied in the theoretical viewpoints. I1,25) Character- 
istic three-peak structures are commonly observed in the 
mid-infrared region around 0.1-0.4eV, where the highest and 
second-highest peaks are called the dimer and Hubbard exci- 
tations, respectively. The lowest peak is sensitive to the band 
width, and is not assigned yet. As shown in Fig. 2, the opti- 
cal absorption peaks in the present theory are classified into 
the inter-MO excitations (peaks C and D), the Hubbard exci- 
tations (peaks B and E), and the intra-dimer excitation con- 
nected to the collective excitation (peak A). In the case where 
the DM and CO phases compete with each other, as examined 
in the present paper and expected in k-(ET) 2 Cu 2 (CN) 3 , ex- 
citation energies for the peaks C and D are higher than that 
for the peak A. A possible interpretation for the experimen- 
tal spectra is that the highest and second-highest peaks in the 
mid-infrared region correspond to the inter-MO excitations 
and the Hubbard excitations, respectively. Strong polarization 
dependence seen in the peak A may rule out the interpreta- 
tion that the lowest peak, which shows a weak polarization 



dependence, is attributed to the intra-dimer excitation. In the 
materials which are far from the DM-CO phase boundary, the 
intra-dimer excitations are hybridyzed strongly with the inter- 
MO ones and may contribute to the spectra around the mid- 
infrared region. 

Recently, Itoh and coworkers find out a new peak in the 
terahertz region around 4meV in k-(ET) 2 Cu 2 (CN)3. 16) This 
peak is only observed in the y polarization. With decreasing 
temperature, the peak intensity increases markedly, in contrast 
to other peaks. This temperature dependence is strongly cor- 
related with the remarkable increasing in the dielectric con- 
stant. 8 ' This newly finding peak in the terahertz region is a 
plausible candidate of the collective charge excitation. It is 
supposed that k-(ET) 2 Cu 2 (CN)3 is located near the phase 
boundary between the DM and polar-CO phases, and below 
25K where the dielectric constant takes a peak, inhomoge- 
neous small polar-CO domains grow up. The increasing of 
the terahertz -peak intensity is attributed to development of 
the polar-CO region in low temperatures. We further expect a 
coupling between the collective excitation and spins through 
the S ■ SQ x Q x -type interaction in "Hj, 9) which gives rise to 
the competitive relation between the polar CO and spin cor- 
relation. When the observed terahertz peak is attributed to the 
collective excitation in the polar CO region, it is supposed that 
increasing of the spin fluctuation shifts the system to the phase 
boundary and reduces the collective-mode energy. More de- 
tailed theoretical studies including inhomogeneous effects are 
required to confirm this scenario. The direct observations of 
the collective mode by the electron-energy loss spectroscopy, 
or the inelastic x-ray scattering are helpful to a comprehensive 
understanding of the charge dynamics. 
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